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ABSTRACT 



Context. Knotty structures of Herbig-Haro jets are common phenomena, and knowing the origin of these structures 

is essential for understanding the processes of jet formation. Basically, there are two theoretical approaches: different 

types of instabilities in stationary flow, and velocity variations in the flow. 

Aims. We investigate the structures with different radial velocities in the knots of the HL Tau jet as well as its unusual 

behaviour starting from 20" from the source. Collation of radial velocity data with proper motion measurements of 

emission structures in the jet of HL Tau makes it possible to understand the origin of these structures and decide on 

the mechanism for the formation of the knotty structures in Herbig-Haro flows. 

Methods. We present observations obtained with a 6 m telescope (Russia) using the SCORPIO camera with scanning 

Fabry-Perot interferometer. Two epochs of the observations of the HL/XZ Tau region in Ha emission (2001 and 2007) 

allowed us to measure proper motions for high and low radial velocity structures. 

Results. The structures with low and high radial velocities in the HL Tau jet show the same proper motion. The point 

where the HL Tau jet bents to the north (it coincides with the trailing edge of so-called knot A) is stationary, i.e. does 

not have any perceptible proper motion and is visible in Ha emission only. 

Conclusions. We conclude that the high- and low- velocity structures in the HL Tau jet represent bow-shocks and Mach 

disks in the internal working surfaces of episodic outflows. The bend of the jet and the brightness increase starting 

some distance from the source coincides with the observed stationary deflecting shock. The increase of relative surface 

brightness of bow-shocks could be the result of the abrupt change of the physical conditions of the ambient medium as 

well as the interaction of a highly coUimated flow and the side wind from XZ Tau. 

Key words. ISM - Stars: winds, outflows - ISM: jets and outflows - ISM: individual objects: HH 151 



1. Introduction 



The Herbig-Haro jet from HL Tau (HH 151) was first 
noted bv lMundt &: Fried! (|1983[ ) and became one of the first 
discovered jets from young stellar objects (YSOs). The HH 
flow associated with HL Tau has an exciting history and 
even today the origin of some its features is not clear. For 

S^ I example, it starts as the relatively faint and narrow jet 
$H ■ from HL Tau with a position angle of about 51° but at 

_^.' a distance of about 20" from the star it abruptly changes 
direction by about 14° and becomes much more brighter 
and wider. This strange behaviour was the main reason for 
several authors to consider this flow as an intersection of 
two jets: one from HL Tau and another from the hypothet- 
ical source VLA 1 (Mundt ct aL 1987). But sub- millimeter 
(JMoriartv-Schieven et al.ll2006( ) and various infrared obser- 
vations of this region did not confirm the existence of any 



* Based on observations obtained with the 6-m telescope of 
the Special Astrophysical Observatory of the Russian Academy 
of Sciences (SAO RAS). The observations were carried out 
with the financial support of the Ministry of Education and 
Science of Russian Federation (contracts no. 16.518.11.7073 and 
16.552.11.7028) 



source in the position of the VLA 1. Presently the HH 151 
jet is considered as a single fiow from HL Tau itsclL In this 
case one should assume that the abrupt changes in its mor- 
phology and direction can be produced by variations of the 
conditions in the surrounding medium. 

On the other hand, the neighbouring pre-main-sequence 
(PMS) star XZ Tau is surrounded by an elliptical shell of 
dimensions 1.5 ^x2', re vealed by the high-resolution imaging 
in 13CO(1-0) (IWelch ct al. 2000). Portions of the wall of 
this bubble can be seen als o in the optical. B esides, XZ Tau 
drives a wide HH outflow (|Krist et al.lll999t ). The relation 
between this flow and the bubble around XZ Tau is not 
clear yet, but we see that HL Tau is situated at the edge of 
an expanding shell of gas and dust, probably driven by the 
quasi-spherical wind from XZ Tau. 

Another interesting feature of the HL Tau jet is that 
its knotty structure (typical for nearly all jets) also be- 
comes prominent only after the bending and brightening 
described above. The Fabry-Perot (FP) interferometry of 
these knots revealed th eir complex kinematical structure 
(jMovsessian et al]l200'<1) . with low and high radial velocity 
components of distinctly different morphology. In addition, 
these knots also demonstrate perceptible p r oper motions, 
which were first studied by iMundt et al.l (|1990[ 1. Other 
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dat a, presented in the pa pers of iMovsessian et aLl (|2007[ ) 
and lAnglada et all (|2007t ). showed a certain difference be- 
tween the proper motions (PM) of the knots, obtained in 
Ha and [S ll] emission hnes. 

There are various theoretical approaches to ex- 
plain the formation of knotty structures in the jets 
(jReipurth fc Ballvl l200lt) . Usually either various types 
of instabilities in the j ets - Kelvin-Helmh oltz instabili- 
ties in hvdrodvn amical (iMicorio et alj 119981) or magneto- 
hvdrodvna mical (lOuved et alJ 200 j) models, thermal in- 
stabilities ([Gouveia dal Pino fc Qp hcr 1990) or time vari- 
ations in the vel ocity of the flow that ge nerate internal 
working surfaces (jRaea &: Noriega-CrespoJ[l998.) . are con- 
sidered. To prefer one of these approaches, we need more re- 
fined observational data. All previous proper motions (PM) 
investigations of HH jets were performed in narrow-band 
images and represent the motions integrated through the 
full profile of the emission line. We attempted for the first 
time to measure the PM in the HL Tau jet for structures 
with different radial velocities and distinct morphologies, 
taking into account the possibility of FP scanning inter- 
ferometry to obtain the velocity-channel images with very 
high spectral resolution. 

Thus, in this paper we present proper motions of the 
structures in the HL Tau jet discussed in our previous paper 
(jMovsessian et al]|2007[ ). obtained from scanning FP obser- 
vations of this region in 2001 and 2007. We also discuss a 
scenario in which the formation of tilted bow-shocks and 
bending of the jet could be the result of interaction of the 
HL Tau jet and the XZ Tau wide flow. T his possibility was 
also p roposed in the our previous paper (|Movsessian et al.l 



also p r 
[2001 . 



2. Observations and data reduction 

Scanning FP observations in the Ha emission line were 
carried out at the prime focus of the 6 m telescope 
(SAO, Russia) on 10 February 2007 (seeing was 1 - 1.5"). 
We used a Queensgate ET-50 etalon, operating in the 
501st order of interference (at the Ha wavelength), placed 
in the coUimated beam of the SCORPIO focal reducer 
(|Afanasiev fc MoiseevI [2005[) . The SCOR PIO cap a bilitie s 
for FP interferometry were described by I MoiseevI (|2002D . 
An interference filter with FWHM = 15A centred on the 
Ha line cut the desired spectral range. The detector was 
a 2Kx2K CCD EEV 42-40 operated with 2x2 pixel bin- 
ning to reduce the read-out time, so 1024 x 1024 pixel im- 
ages were obtained in each spectral channel. The field of 
view of about 6' was sampled with a scale 0.36" per pixel. 
The spectral resolution was estimated as FWHM « O.SA (or 
PS 36 km s^^) for a range of AA= 13A (or sa590kms^^) free 
from overlapping orders. The number of spectral channels 
was 36 and the size of a single channel was AA «0.36A 
(«16kms-i). 

We reduced the interferometric data with an IDL-based 
software developed at SAO (Moiscev 2002). After the pri- 
mary data reduction, the subtraction of night-sky lines, and 
wavelength calibration, the observational material consists 
of "data cubes" in which each point in the 1024 x 1024 pixel 
field contains a 36-channel spectrum. 

Proper motions were determined by comparing the ve- 
locity channels, corresponding to —50 km s~^ and —150 km 
s""'^, obtained in 2001 with the same telescope and optics 



but another detector (see IMovsessian et al.l 120071 ) . and the 
new ones of 2007 (all radial velocities used in this paper 
are heliocentric). To increase the signal/noise (S/N) ratio 
we summed the three channels in each cube around the 
above mentioned central velocities. The velocity range in 
both cases was equal to 31.5 kms ~^. Then we determined 
the PM for both high and low radial velocity structures. 

The first step of a PM measurement from FP data is to 
reduce the data cubes of both epochs to a common reference 
system. For the HL Tau field this procedure has certain dif- 
ficulties: the reference stars are scarce, and the greater part 
of these stars are embedded in reflection nebulae. Indeed, 
we were forced to use nine reference stars, including such 
nebulous objects as HL Tau, XZ Tau and LkHa 358, which 
definitely diminished the precision of astrometry. The rms 
of the differences in the positions of reference stars is about 
0.24". 

The "diffe r ence square d" method, d escrib ed in 
iHartiean et"an (|200lD and iReipurth et all (|2002t) . was 
used. This algorithm is well- adapted for the measurements 
of extended objects and works in a selected rectangular 
region in the images of both epochs, finding the shift in 
the plane of the sky, which minimizes the sum of squared 
differences between the intensities in the defined region. 
For comparison the proper motions were also determined 
using the barycenters of the knots. For the HL Tau jet 
the systematic differences between both methods were 
negligible for the compact high-velocity knots, but for 
the low- velocity bow shocks they became significant. 
Therefore, only results from the difference square method 
were used for the proper motion estimates for the bow- 
shocks. We estimate the total uncertainty of the position 
of the correlation peak and astrometry to be about 0.3" 
for the compact knots and about 0.4" for the extended 
ones. 



3. Results 

As was already des cribed in our previous paper 
([Movsessian et al.l l2007l ). the HL Tau jet consists of the 
high-velocity narrow jet with knotty morphology and the 
low- velocity bows, which are visible only in Ha and cross 
the jet, being located in front of the high-velocity knots. 
They exhibit the greater transverse widths at the greater 
distances from HL Tau and are not axisymmetric, being 
elongated in the SE d irection. Their m orphology is best 
visible in HST images ([Krist et al.ll2008D . 

Th e knots in HL Tau jet were described bv lMundt et al.l 
([1987D and are labelled from A to E. The connected bow- 
like structures were detected for knots B, C, and D only. For 
knot A no visible variations of the morphology in various 
velocities were detected; the profile of Ha emission in this 
knot is wide but shows no split into two velocity compo- 
nents. Knot E is visible only in the channels corresponding 
to high radial velocity. Thus, we estimated the PM of low 
and high radial velocity structures only in the knots B, C, 
and D. However, for the sake of completeness we measured 
the PM for knot E in the high-velocity range and for knot 
A using the full line profile. 

The results of our PM measurements are presented in 
Fig.[l] where the two images of the HL Tau jet, correspond- 
ing to low and high radial velocities, with the vectors of PM 
of the knots, are shown. The numerical results, again sepa- 
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Fig. 1. Proper motions of the structures in the A, B, C, D, and E knots of HL Tau jet, corresponding to high (a) and 
low (b) radial velocities, are shown by vectors plotted in the images in Ha emission corresponding to —150 kms~^ and 
—50 kms ~^, respectively. The scale of the vectors is indicated by the arrow at the top of the left panel. Rectangles show 
the areas used for the difference square method. 



Table 1. Proper motions in the jet of HL Tau 



knot AQf(arcsec) 



A(5(arcsec) 



Vt (km s-^) P.A. (dog) 







full velocity range 






A 


13.6 


10.9 
high radial velocity structures 


60 ± 40 


45 ± 19 


B 


16.4 


14.7 


162 ± 40 


47 ± 12 


C 


20.3 


18.6 


157 ± 40 


46 ± 15 


D 


24.9 


27.7 


155 ± 40 


47 ± 19 


E 


34.7 


34.6 
low radial velocity structures 


150 ± 40 


47 ± 19 


B 


17.1 


13.7 


154 ± 50 


44 ± 15 


C 


19.8 


17.1 


167 ± 50 


35 ± 24 


D 


27.6 


26.6 


161 ± 50 


36 ± 17 



Column 1; Knot (nomenclature from lMundt et aU (| 19901 )): Columns 2 and 3: Right ascension and declination offsets from the 
position of HL Tau; Column 4: Proper motion velocity, assuming a distance of 140 pc; Column 5: Position angle. 



rately for the low and high radial velocity structures (except 
for knot A) , are given in the Table [TJ 

The most obvious and interesting result is that the 
structural components with low and high radial velocity 
both have very similar values of proper motions. We will 
discuss this in more detail below. R is also worth men- 
tioning that these new data confirm the tendency of the 
PM vectors of low-velocity Ha features to turn increas- 
ingly north with increasing distance from HL Tau, found 
bylMovscssian ct al. (2007). 

When considering the motions of the knots between two 
epochs their shift is obvious except for knot A. Rs behaviour 
can be described as the progressive elongation in the direc- 
tion of its motion, leaving its trailing edge nearly fixed. 
This is clearly seen in Fig. [^ where the images of knots A 



and B in the two epochs, integrated by full line profile, are 
presented. 

For further clarification we show in Fig. [3] the intensity 
slice along the HL Tau jet in integral Ha e mission. The 
nomenclature of the knots is the same as in iMundt et al.l 
(|l990l ). The shift of the knots between two epoches is ap- 
parent except at the beginning of knot A where the abrupt 
brightening of the flow begins. This point can be considered 
as stationary relative to HL Tau. The change of the relative 
brightness of knots A and B between 2001 and 2007 should 
also be noted. 

Comparing the present measurements with the pre- 
vious ones, we see that the directions are quite similar 
throughout the literature, though the values differ signif- 
icantly. The present PM values are more similar to the re- 
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Fig. 2. Appearance of knots A and B in the HL Tau jet 
in 2001 (contours) and 2007 (greyscale). The progressive 
elongation of knot A and the obvious shift of knot B can 
be compared. 




arcsec 

Fig. 3. Cuts along the HL Tau jet for the two epochs of 
observations (full line: 2001, dashed line: 2007), with the 
designations of the knots shown below. The intensities of 
both cuts are normalized by knots C and D. Shifts of the 
knots are apparent except for the beginning of knot A. 



su its of Anglada et al.l (|2007) than to the hig h er est imates 
of iMundt et al.l (|l990t ) and iMovsessian et al.l (|2007D . Very 
probably this discrepancy results from the unavoidable sys- 
tematic errors of various origin, which were briefly discussed 
also in SectlJl In any case, this inconsistency cannot affect 
the basic results and the subsequent discussion in any way. 



4. Discussion 

4.1. Multiple working surfaces in the HL Tau jet 

The complex kinematical picture of knots in the HL Tau 
jet speaks about the existence of two distinct shock struc- 
tures in each knot. These structures differ not only by 
their radial velocity but also by excitation, because bow- 
shape structures are visible only in Ha e mission and high- 
veloc i ty knots are much brighter in [S ii] (Movs essian et all 
120071: iKrist et al.ll2008[ ). This picture strongly resembles 
separation by excitation of bow-shock and Mach disk 
structures in the terminal working surfaces of the jets 



(jRaea fc Norieea-CresDO.! [19981 ) . This division of shocked 
regions was indeed was found through narrow-band im- 
agery (jReipurth & Heathcotc 1992). Moreover, with the 
usage of spectral image ry, these regions w e re also sep- 
arated kinematicallv flHartigan et al.l l2001t iMorse et al.l 
119921 : IMovsessian et al.l l2000t) . But, for the jet knots even 
the doubtless splitting of two shocked structures - both by 
radial velocity and by excitation - leaves some uncertainty 
in choosing between the theoretical mechanisms of the knot 
formation. Here, in our opinion, the proper motion values 
combined with the radial velocities can play the crucial role 
in clarifying this situation. 

As we already showed, the proper motions of the two 
distinct radial velocity structures in each knot of the 
HL Tau jet have the same values. To analyse this picture we 
suggest considering the knots as the internal working sur- 
faces formed in a time-dependent flow, i.e. in the same way 
as the terminal working surfaces. In this case they should 
contain the two principal shock regions: the "reverse shock" 
(or Mach disk) that decelerates the supersonic flow, and the 
"forward shock" (or bow-shock) that accelerates the am- 
bient medium. The pattern speed of these two structures 
should be the same, in contrast to the fluid speed of emit- 
ting particles. 

Indeed, our data allow us to obtain for each knot both 
the fluid speed, which corresponds to the velocity of emit- 
ting particles, and the pattern speed of the visible struc- 
tures. Tangential velocity is the projection of pattern veloc- 
ity onto the plane of the sky: Vt = Y pattern cos a, where a 
is the angle between the outflow axis and the plane of the 
sky. Accordingly, the radial velocity of the knot corresponds 
to the projection of the velocity of emitting particles (fluid 
speed) on the line of sight: V^ = Y fluid sin a. 

The ratio C, of pattern speed to fluid speed was for the 
first t ime determined for HH 34 jet by lEisloffel fc Mundn 
(|l992l ). but the kinematical structures in the knots were not 
separated in their work; only the spectral imagery technique 
allows this. 

To obtain the correct estimates for the jet absolute ve- 
locity and inclination (a angle) we need to use the struc- 
tures where C = 1- Foi' example, in a Mach disk by definition 
the fluid and pattern velocities have nearly the same values, 
so C, should be near unity. Thus, combining the tangential 
and radial velocities of high- velocity compact knots (which 
we assume to represent Mach disks), we find that a — 50° 
± 10° and the mean value of the absolute velocity of the 
jet is about 250 km s^^. 

Previous observ ational estimates for a put its value in 
the r ange 40°-60° (|Anglada et al.l 120071: IMovsessian et al.l 
|2007| ). On the other hand, models of the HL Tau circum- 
stell ar dusty environm ent provide an angle of the flow about 
43° (jMen'shchikov et al. 1999). As one can see, the new val- 
ues agree well with the previous estimates. 

Now we can consider the low radial velocity, bow-shape 
structures, which we assume to be the shock waves that 
entrain particles and rise the temperature of the ambi- 
ent medium. Knowing the inclination angle, we obtain the 
mean value ^=2.6 for the three bows in the HL Tau jet, 
which strongly diffe rs from the comp act knots. 

As was shown bv lHartiganI ()1989t ). the ratio of radial ve- 
locities in shocked structures depends on the flow and am- 
bient medium densities as well as on the jet velocity itself. 
From the same model, it follows that the relative surface 
brightness of a Mach disk and a bow-shock also depends on 
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the ratio of the jet and ambient medium densities (n-,- and 
lia correspondingly) , and brighter bow-shocks are formed in 
denser jets. For knots within the HL Tau jet, the relative 
surface brightness of bow-shock and Mach disk is nearly 
equal, so the ratio na/n^ is a little more than unity. Thus, 
the sudden increase of brightness of the jet and appearance 
of bright bow-shocks, starting from the 20" distance from 
the source, could be understood as the result of an abrupt 
change of the physical conditions in the ambient medium, 
when the flow reaches the region of lower density formed 
by the poorly coUimated outflow from XZ Tau. 

Of course, one cannot exclude the possibility that this 
brightening can be the result of the lower line-of-sight ex- 
tinction after the deflection point. Ho wever, both in our 
images and in Fig. 4 from the paper of iKrist et al.l (|2008[ ) 
we can see that the narrow jet is visible from the star and 
almost to the deflection point; no Ha emission low- velocity 
structures exist in this range. 

4.2. Jet - wind interaction? 

Now we can consider the unusual behaviour of the HL Tau 
jet, gathering all its unusual features. Some of them were a l- 
ready discussed in our first paper (Movses sian et al.ll2007l) . 
In summary, we have the following points: 

• The bend of jet at the distance of about 20" from the 
source star; 

• the abrupt increase in the brightness of the jet, starting 
from the same position; 

• position of this point is not shifted between the two 
epochs of observation; 

• appearance of low radial velocity structures only after 
the bending of the jet; 

• these structures have the shape of non-axisymmetric 
bows with brighter SE edges; 

• transverse velocity gradient is observed in these low- 
velocity components. 

We suppose that these phenomena are mutually inter- 
related and, as was mentioned above, that they can result 
from the abrupt change in the physical condition. 

Deflection of the HL Tau jet to about 12° could oc- 
cur from the glancing collision between the jet and the 
cloud m aterial, located ap proximately at 20" distance from 
HL Tau (|Krist et al.ll2008D . as in case of HH 110, where the 
jet collides with a dense clump (Reipurth et al. 1996). But, 
as was mentioned above, the relative brightening of bow- 
shocks (which is observed in this case) takes place when the 
flow reaches the region with lower ambient density. 

In this connection it is important to consider again knot 
A. Unlike the other knots, its edge, nearest to the star, has 
no perceptible proper motion (Figs. [2] and O; on the other 
hand, the knot is detectable only in Ha emission (Fig. 2]). 
We are inclined to believe that the tail of knot A represents 
a so-called deflection shock, formed at the dense wall cre- 
ated by the wind from XZ Tau. Both its stationarity and 
morphology have a great similarity to the f eature that was 
recent ly discovered in the HH 46/47 jet bv lHartigan et al.l 
(1201 ll) fcompare the ir Fig. 14 and the Fig. 4 from the paper 
of lKrist et al.l ([2003)) and considered by them as the same 
kind of deflection shock. 

On the other hand, the observed bending of the jet can 
be produced by the ram pressure of a side wind, as was 




Fig. 4. Superposition of the HL/XZ Tau FP images in the 
monochromatic [S ii] A6716A line emission (greyscale) and 
low- velocity Ha emission (contours). The Ha observations 
were performed in Sept 2001 (see Sect.[2l, the [S ii] data 
were taken in Jan 2003 (jMovsessian et aLll2007l ). Knot A 
(shown by arrow) is visible only in Ha and appears as a 
continuation of the Ha emission arcuate structure west of 
XZ Tau. 



shown in ICiardi et al.l (I2008D. An expa nding shell around 
XZ Tau, detected bv lWelch et al.l (|2000D . can definitely play 
the role of such a side wind. As a result of this interaction, 
the increase of shock energy release is manifested by the 
brightening of the jet. The formation of elongated wings 
of the bows in the direction of the XZ Tau bubble (to the 
SE) as well as the transverse velocity gradient detected in 
bow-shape structures support this idea. In any case, the 
transverse velocity gradient is too large to be explained by 
jet rotation. 

One cannot exclude the possibility that both mecha- 
nisms can play role in the observed behaviour of the HL Tau 
jet. Using the model described by Can to fc Ragal (|l995l) we 
estimated, h owever, that t he expansion velocity of XZ Tau 
bubble (, Welch et al.l 120001) is insufficient for the observed 
abrupt bend of the HL Tau jet. But, as was already men- 
tioned above, it can be responsible for brightening of the 
jet. 

We also recall that accord ing to our previous measure- 
ments ([Movsessian et al]l2007[ ) , confirmed also by other au- 
thors, in the SW direction from XZ Tau the definite turn in 
the PM vectors of the knots in th e so-called "Ha jet" is ob - 
served. As was also suggested bv iMovsessian et all ()2007[ ). 
this can be caused by the interaction between this flow and 
the south-western side of XZ Tau bubble. These two cases 
of the unusual behaviour of the jets around the XZ Tau, 
in our opinion, argue that a real interaction is taking place 
between the jets and the expanding bubble of this star. 

5. Conclusion 

For the first time the proper motions of the knots in the col- 
limatcd jet were measured with the scanning Fabry-Perot 
interferometer. This allows us to make certain suggestions 
and conclusions. In the end we summarize our findings in 
the following way: 
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(1) Our PM measurements show that the two structures 
distinguishable by radial velocity that exist in the HL Tau 
jet have the same tangential velocities, with a typical value 
of the order of 160 km s~^. 

(2) The angle of the flow to the sky plane, determined 
from the radial and tangential velocities of high-velocity 
structures, is 50° and the total space velocity of the flow is 
about 250 km s~^. 

(3) We are inclined to assume that the different velocity 
structures in each knot of the HL Tau jet are detecting two 
principal shocks, namely bow-shock and Mach disk. 

(4) Taking into account all observed features, we con- 
clude that the knots in the jet originate from the episodic 
velocity variations in the HL Tau outflow, and each of them 
forms internal working surfaces with two principal shock re- 
gions. 

(5) The strange and unusual behaviour of the HL Tau 
jet, starting from ~20" from the source, is connected with 
the abrupt change of physical conditions in the ambient 
medium. Namely, the jet enters the zone where the matter 
is blown out by the wide wind from XZ Tau, creating a 
deflection shock on its border where the swept-out matter 
formed the dense shell. 

Recent narrow-band observations of the knots in HH 34 
bvlHartigan et ahl ()2011t) demonstrated that those can also 
be resolved to the high- and low-excitation regions, and 
that in each knot both shock regions have the same proper 
motions. This interesting result can be considered to sup- 
port our interpretation of the high and low radial velocity 
structures in the HL Tau jet. 

We hope that this work, together with other new results, 
will give an impulse to develop new theoretical models of 
the collimated outflows, which will take into consideration 
all modern observational data. 
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